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The relative retardafcion between the components of a stream of 
light produced by the passage through a crystalHne plate can, as is 
well known, only be determined in finite tei-ms in a limited number 
of special cases. Tn general it is necessary to be content with an 
approximate solution, and those hitherto published have, as far as I 
have ascertained, never been carried beyond terms of the second 
order with respect to the sine of the angle of incidence of the light, 
while they do not readily lend themselves to a further approxima- 
tion. The following method of dealing with the problem, which, 
except for the labour of calculation, can easily be extended to 
terms of any order, may then be of practical use, and, as leading to 
an interesting relation between the corresponding terms in the 
development of the roots of a certain biquadratic equation, of some- 
interest. 

1, Let the pla.ne of the paper represent the plane of incidence, OT' 
represent the normal to the front of the incident plane wave, OP and. 
OP' the normals to the fronts of the two corresponding refracted 
waves ; then if OM, OF, and OF' represent the spaces these wa^es 
would traverse in unit time, planes through M, F, and F' perpendi- 
cular to OT, OP, and OP' respectively will i-epresent the positions 
that the fronts of the waves would occupy in. unit time after leaving 
O, and these planes will, by Huygens' principle, intersect the 
surface of the crystalline plate in one straight line, projected on the 
plane of the figure in the point 0. 

Let OT, OP, and OP' meet the second surface of the plate in the' 
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points Tj P, and P' respectively, and through P and P' draw PE, 
P'E' perpendicular to OT, meeting it in the points E and E' ; through 
E' draw E'T' parallel to the surface of the plate, meeting PE in T'. 

Draw OD normal to the plate, meeting the second surface in D, 
and let the wave fronts ON and OW meet this normal in the points 
Q and Q'. 

Then the ti^angles OPD, COE", CQO are similar, as are also the 
triangles OP'D, CON', CQ'O, and the triangles T'E'E, COM. 

The relative retardation, measured in time, of the two waves after 
both have traversed the plate, is represented by 
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Hence if the axis of Z be normal to the plate, and 

ZX-hmY+%Z = 1 ZX-hwY + ^Z = 1 

are the equations of the refracted waves in unit time after passing 
through 0, their relative retardation, measured in length in air, after 
both have traversed the plate is 

A = ^T(% — %), 

where v is the propagational speed in air, and T is the thickness of 
the plate. 

Also i being the angle of incidence, and w the azimuth of the 
plane of incidence with respect to that of XZ, 

I = sin % cos wivy m = sin i sin wlv, 

2. In applying this result to the case of a plate cut in any manner 
from a biaxal crystal, let the surface of the plate on which the light 
is incident be taken as the plane of XT, the plate lying on the side of 
Z positive, and let the plane of XZ be taken so as to contain the 
least axis of elasticity of the crystal. 

Let Ox, 0?/, 0^ be the axes of elasticity, and the angle yoY = 0, 
the angle zo7i = y^, then the transformation from the axes of elasti- 
city to the new axes may be effected by the following successive 
transformations, each in one plane : — 

(1) Through an angle 0, in the plane of xy^ from Oo?, Oy to 

0^1, 0^. 
('2) Through an angle x, in the plane zxi, from 0^, Oa^i to OZ, OX. 
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The formulge for these transformations are — 

m = Xi cos 0-- Y sin 0, y =^ sci sin 04- Y cos 0, 

a?i = X cos x+Z sin x? ^ = —X sin x+Z cos x> 

from which we obtain 

a? = X cos cos X— Y sin 0+ Z cos sin x, 

2/ = X sin cos x+ ^ cos 0+Z sin sin x? 

^ = — X sin x+Z cos X- 

ITow the equation to the wave surface referred to the axes of 
elasticity is 



• „2 7.2 I -.2 „2 ~~" ^5 



r^ — o,^ r^ — 6^ r^ — c^ 

a, 5, c "being the principal wave velocities, and the condition that the 

plane 

Ix+my-j-nz = 1 

should be a tangent plane to it, is obtained by eliminating p between 
the equations 

-r V ' » / a^—p^ Ir—p^ cT'-'p^ 

Hence the condition that in the new system of co-ordinates the 
plane 

ZX+mY+^Z = 1 

should touch the wave surface is found by eliminating p between the 
equations 

p^ = if? -j- m^ 4- 5*^)""', 

(I cos cos x—m sin 0+^^ cos sin x)^ 

(f sin cos X + m cos 4- *?* sin sin x)^ Q sinx-~"^cosx) ^ _ ..x 

A ^ L • - 2 — 0« • (1). 

f 5>5 ,p^ C^—p^ ^ \ 

The result of this elimination is a biquadratic in n^ which, from 
the nature of the problem, has two real positive and two i^eal nega- 
tive roots, and if ni, n^ are the positive roots of the biquadratic, the 
relative retardation required is 

A = -yT (%— W2). 

3. Before proceeding with the general case, these results may be 

applied to certain simple cases : — 

n 2 
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(1) Let the plate be cut from a uniaxal crystal, then writing^ 
?> = a, equation (i) gives the two equations 

jp^ = a^ and 

(p^— c^){(Zcosx-f wsin%)^4-m^} + (jp^— a^) (I sinx— ^cosx)^ = 0, 
and the values of n are given by 

a'^(P-\-m'^+n^) = 1, 

c^(P+m^-\-n^)—l-\-(a^'-c^) (I sin x— '^^ cos x)^ == 0, 
whence 



• o • 



1 / „ sm'' ^ 
% = -A / l — d — 5-j 
a V V 



Ui = { V'(a^cos^x + ^^sin^x) (l — d^smHlv^)—c^(a^'-'G^) sin^x gos^w sin^'/^^ 

+ (a^— -r) sin X cos X cos w sin i'/v} 4- (a^ cos^x+ ^^ sin^x) 

and 

A ^/^~ — a'^ sin^ ^ (a^ — c^) sin x cos x cos w s In ^ 
T ~ a a^ cos'^ X + c^ sin*-^ X 



y(a^ cos^ % + <^^ sin^ x) (f^^—c^ sin^ ^) — c''^(a^ — c^) sin^ x cos^ zv sin^ ?-' 



2 „-;^2 



crcos^X + ^^i^ X 

(2) Let the plate be cut from a biaxal crystal perpendicularly to 
the mean line ; then taking the axes of elasticity as the co-ordinate 
axes, the biquadratic in 7i becomes 

- {(P + c^) P + (c' + a') m^ 4- (a^ + ¥)n^} + l = 0, 
Cir 

and ^1, ^2 being the positive roots of this equation 

(^-^2)2^252 ^ (^2^^2)_52 (^^2^^2) ^2_^2 (2,2^^2) ^^2 

and 

aW-fT^ = (^a^-i-&^) i;'"— {6^(c^ + tiO cos^i(;+<^^(&^4-c^) sin- 21;} sin^^ 



— 2a6 V (t?^ — c*"^ sin''^ i) { -y^ — ( 6^ cos^ w + <^^ sin^ tc') sin^ i } 
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4. Returning to the general equation (i), write 

A' = cos sin X? A = cos cos x cos ta— sin <j6 sin w^ 

B' =: sin sin x, B = sin <p cos % cos w + cos sin w, 

C = cosx) = —sin X cos w. 

Then the equations between which p has to be eliminated become 

1 „ sin i 



and since 

the result of the elimination becomes 
A'V+2AA' ^n+A' —~] (hVn'+hV ~~- V-c' 

+ ( B' V + 2BB' — - 7^ + B'^ —l^] (c'aV + &o? ~™ --c'-o? 



• # 



;,'.2«'\ / «-;^.2/ 



+(C'V+2C0'^-w+C^— ~ (a^6V+a^&^— T^-a^-P) + l = 0; 



V v^ J \ v^ 



and multiplying out and arranging the terms, this equation may be 
written 

sin^' „ . / sin^n „ L sin^ . \, sin^ A sin^^ 



. _ sin^ „ , / sm~^\ „ /, sin^ , ^ sin"*\ 
?i*--2ai n^-\-\ C0 + C2 — T- 1^—2 61- -I- 63 — T-h^ + <^o + t^2 



where 



sin^^ ^ .... 

4-0,4 — ^~ =0 (11), 

- ^ _ ^AX^' + cQ _ SA^6V 

^' ~ 2 A'^6^c^ ' "'' ~ 2 A'^6V ' ""' ■" 2A'^6V ' 

^^ ~ ~~2A¥?" ' ^3 - ^1 - --^^.-^ ' 

_ 2A'X&'+o^) _ 

Co — 2A^V~ ' ^^ "~ '^■^'^*' 
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Now the roots of equation (ii), p, q, r, s will be functions of sin i/v, 
and expanding these in series proceeding according to powers of this 
quantity, we may write generally 

sin^ , sin^^' , 

sin i , sin*^ i , 
q = ^0 + 2^1 1-^2 — T-+ •. •• 

sin i , sin^ i , 
r = ro-fn \-r2—:r~+ •• •• 

V V" 



sm^ sin^^ 

V v 



But 



• • 



sm % „ . sm'^ i „ , sin i . , sin^ i 



2p = 2di , 2pg = C0 + C2 — 7- , 2pgr = 261 \-2bz — -— , 



sin^ i . sin^ ^ 



Spgrs = ao-\-a2 — y- + ai 



v" v^ 



and therefore equating the coefficients of like powers of sin^'/t? on the 
two sides of these equations, the coefficients p, q, r, 5 may in general 
be determined in succession by means of linear equations, provided 
we can determine Pq, q^, Tq, Sq, which may be at once done, since they 
are the roots of the equation 

n^ -\- Con"^ -{- CLo = 0. 

Moreover, when this method is practicable, the roots of the 
biquadratic take the following form 



where 



sin'^^' sin^^ 

TTi — 0-0 + 5^2 2 tOCi 7 h •• • . , 

V V 

TTg — po-f-Pz ^ \-pi 7 h • • • • ? 

V V 

sin i sin^ i sin^ i 
Pi — 71 r 73 — s — r 75 — r- + • * • • > 

V V V 

sin i sin^ i sin^ i 
p2 = h 73 — o 73 — I — !-..•. 

For suppose that this is true as far as the terms involving 
sin'*""^^7^^~\ and let a, ^, 7, ^ denote the sum of the terms in ttj, TTg, 
/?i, p2 respectively of an order less than n, then we may write 
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and neglecting tlie products of tlie /t's witli one anothei', and with ^ 
and ^j since sucTi products will introduce terms that we shall not 
require, we have 

Hence if p», g«j r,j, % denote the coefficients of shiHIv^ ^^Pi f> **j ^ 
respectively, and [ ]» denote the coefficient of bihHIv^ in the expres- 
sion within tbe vinculum, the equations that determine these qnan« 
titles are — 

I) If » = 2m 

miizm—pzm) + fioiszm'-nm) = C2,»+ [^^ + ^^— («^/ + a)^— 27a]2»i 

== Itggi, say, 

1 

= Sgjje, say 5 
and from the first and third of these equations, we have, unhs$ 

which relations hold whatever m may be. 
(2) li n=- 2ot+1> 

^0 (3^1 ?» + 1 •~*i?2»+l ) 4" ft ( S2m + 1 — ~ **2m+i ) = 0, 

— fi(?(pm+l'i-q2m+l)--0^(?(nm+i+B2m+l) — 262M+i + 2[ai^a 

+ fi^^ — ^B (>Y + ^) ]2W+l 

^^= J- 2m -fly say, 

and from the second and fourth of these equations 

^2w+i = Pmj^h HviJ^i = Tm+u 'Unless a^ = ft^^ 
and when m>0, (I2W+1 = 0, so tliat nm^i == p%mJr\ («^>0). 
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Leaving tlien for the present the case in which a.o' = y8o\ the roots 
have the form given above, and 

OCq — po ^0 — po oco — po 

<ajo — Po , <3£o — po 

whence 

A 

— (Zo—poH- 5 ;7i • — - 2' 



T3 sin^ i j^ R4 + S4 sin^ i 
oto—Po' -0^ «o—/3o ^^ 

the terms of which may be readilj calculated in succession. 

5. It is now necessary to consider the case reserved above, in 
which cco^ = pQ^. When this occurs, the two refracted waves, corre- 
sponding to normal incidence of the light, traverse the plate with the 
same velocity, and hence the plate is perpendicular to an optic axis. 
The case, then, of a plate cut perpendicularly to an optic axis 
requires a special investigation. 

Now for a plate so cut— - 



= 0, cosx = y(&^— c^)/^/(a'^— c^), sinx= ^/{o?—¥')l^{a?—c^), 

^and the biquadratic (ii) has the form 

, -.sin-i o , / r , ,sin^4\ „ /, ,sin^ . , ^sin^i 
m^—2di ' #4- Co +C2 —IT" n?—2 hi ■ h&a — t" I^ 



• o • 



, ,sin"'^ ,sm*^ 
+ ao +a2 — ;r~ + «4 — r" = 0? 



* 4. 

r 



^~. ^ 



where 



di = 63' = cos w • T5 \/(a^ — b^) (¥ — c^) , 



hi = — cosey * v4\/(a^— 6"'')(6^-~-c^), 



2 , , a^c^ , 1 



^/ = —72? C2' = l + a/ = l + "7j- + T4(»^-~&0 (&^ — C^) COS^'M?, 



5^ 



^^0 — T~4> (^% — Tr~~ ' 
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Sln^ 



whence writing n-\-\di —— for n, the equation becomes 
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where 



0, 



Co 



2 1 

T-2, C2 = fi {M + aV— l(ct^--P) (5^ — c'-) cos^iy} , 



\ • — r g — ^ {a^—P) (¥—c^) . cos w, 



<^0 — -^4 
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Let US, in the first place, neglect the coefficient of n, then we 
have — 



o^' 



j^ J sin^^ / , 

2 I — Co — C2— ^-± ^ Co • 



■4^0+2 (C0C2— 2^2)— ::r" + (^2^- 
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4hXa'~h''){h''-c') 



, '6-^ -f a'-^c'^ o 



6^ 



and the roots are ± (t^+z^)? ± (^~/^)) where 



sin^ ^ 



TT = j:?o + p2 — r 4-i?4 



sin'^ i 



V 



V 



" 1 ■ 9 • a *, 



• • 



sm^ . sui^^ 



p = p^ ■ i-p3 — g-+ . , . „ 



V 



sin^ i'^ 



V 



and writing for shortness 



P = ^/(a' 



W) (63-cO, Q = :i (a^^~60 Qf^-o'') cos^t(?-i(6* + aV), 






(a'^-6O(6'--c'0 



1(6^ + aV) cos^ w, 



we have 
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^» = X' -P' = S' ^^ =i (Q-iP'-). i'3 = ^^(PR-PQ+iP). 

Suppose now that the actual roots of the biquadrafcic are — 

then 
(7r4-/>) (^— ^)H-(7r-~/5) (a-7) 4-^/3 + 7^+ (^ + /3) (7+^) = 0, 



sin^i 



^?=^ 



a/37^H- (TT-f-^) 7^ (^ — (%) -\- (jr — fi) ocp (^ — 7) 

-(7r~py{a{3+(7r-^p) (/3-r/)} = 0, 
which, on reducing and simplifying, become 

4:r(a+7)+4/3(a-i-a)+a:H/^^ + 7H^^ = 0, 



sin^^' 



i-7r(a: + 7) (^^ + 5:7) +/)(<% + a) (^7 + 5i^) +afir^B = 0. 

Let 

sin^ ^ sin^ 1 sin* ^ 

a = 5^2 — —-T^'d — 3-4-^4 — 4-+«».« ? 
i;-' 1'^ -y* ^ 

with, similar expressions for /3, 7, S ; then the terms involving sin^ijv'^ 
give— 

<»2 + A + 72+^2 = 0, ^24-72=0? 

. • . 0^2 = --^2 and ^2 = — Aj 
and fi^om the terms involving sin^/'/'^"^? ^'^ obtain — ■ 

0^34-/534-734-^3 = 0, 

i>o(a3 4-73) 4-^1(^24- ^2) = 0, 

CC'i — ft'l — 0, 
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wlience 

^2 = p2 = —72 = —^2 = hlQ^PoJpi) = "p • h, 

and 73 = -^3, ^3 = -A. 

Again, from tlie terms involving sin.^i/v^-"- 

^4 4- A + 74+ ^4 = 0, 
Po(^4 4-^/4) +i>i((%3 + a3) +fli/ = 0, 

as+ft = 0, 

whence ^3 = -ft = -73 = ^a = - ^^ = --— • W, 

and oci = -— fyi, A = — b^. 

Similarly tlie terms involving sin- ijv^ give — 

a5 + AH-75 + 4 = 0» 

j3oPl(a4 + A) + 2 (pipa +1?0P3) ^2 + 2PqOC20C3 = 0, 

i?o(a4—A) +fi(a3— A) +a2^ = 0, or ^4—^84 = 0. 

Hence 

m ^^3 



ai — Pi— ~-fy4 — —04 — = —-— . 534. 

popi jpi 4P 4P^ 

and 0C5 = —75, p5 = ^^6, 

and so on. 

Hence, as far as terms of tlie fourth, order, 

sin^ i sin^ i , sin* i 

a=z — fy = ^2 J- +^3 -~ +Cf4 T~ , 

siii'^i sm^ i , sin^-t 

^ =z —b z= a2 — ^ cx,z— — tai — T- , 






where 



6*-aV __ (h^^a^G^y 



— (&*— aV)^ sin^ w} cos w^ 
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and the difference of the positive roots being 2/3+ a —7, we have, to 
the same degree of approximation, 

sinz sin^^ . , ^ sin^^' . sin^^ 



2pi V2%2 — - — 1-2(^3+^3) — I — V^oci 



tjT V v^ v^ V 



4 



s/(a^—h'^)(¥-~c^) sin ^' , &^~ttV sin^i 

^_- -J^ L^_ ^1.. cOS^t' r- 

, , (a^ — d^y . „ sin^^ , , l^—a?c^ 

+ 8 o ■ ~ ^^ - ^iirio — r. — Vti 



■ 4. • 
X {2{a^—W) (¥-~c^) (&* + aV) cos\d? — (6*-"aV)'^sin^u'}coszt?^-4- • 

6. The proposition on which the aboye investigation depends was 
first suggested to me by an analogous theorem given by McOullagh,* 
in connection with the surface of wave-slowness, f or, as he terms it, 
the surface of refraction or index surface ; in fact, the one may be 
deduced from the other by reciprocating with respect to a sphere of 
unit radius concentric with the surfaces. 

I have since found that ZechJ has employed the same principle 
for the determination of the rings of biaxial crystals, but his method 
of dealing with the biquadratic equation is essentially different from 
that given above, and leads only to the determination of the terms of 
the second order. 

My thanks are also due to Mr. J. L. S. Hatton for some useful 
suggestions that led me to the adoption of the above methods of 
approximation . 

* ' Collected Works,' p. 46. 

t The first pedal of tlie wave-surface is sometimes erroneously called the surface 
of wave-slowness j but, as Sir William Hamilton calls the inverse of the wave- 
velocity the wave-slowness, the inverse of this surface, or the polar reciprocal of 
the wave-surface, is projDcrly the surface of wave -slowness. That this was the 
name given to the polar reciprocal of the wave-surface by Sir William Hamilton 
appears from Lloyd's " Eeport on Physical Optics " (' Collected Works,' p. 122), and 
from McOullagh (' Collected Works,' p. 96), though in his papers he calls it the 
surface of components of normal slowness. 

J ' Pogg. Ann.,' vol. 97, p. 129 ; vol. 102, p. 354. 



